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Abstract 

The thermodynamics of massless (/)^-theory is studied within screened pertur- 
bation theory (SPT). In this method the perturbative expansion is reorganized 
by adding and subtracting a mass term in the Lagrangian. We analytically cal- 
culate the pressure and entropy to three-loop order and the screening mass to 
two-loop order, expanding in powers of m/T. The truncated m/T-expansion 
results are compared with numerical SPT results for the pressure, entropy 
and screening mass which are accurate to all orders in m/T. It is shown that 
the m/T-expansion converges quickly and provides an accurate description of 
the thermodynamic functions for large values of the coupling constant. 
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I. INTRODUCTION 



The behavior of finite temperature field theory at intermediate to large coupling is of 
particular interest due to the upcoming heavy-ion collision experiments at RHIC and LHC 
For years the hope has been that due to the asymptotic freedom of QCD that weak-coupling 
expansion calculations within finite temperature field theory would suffice to describe the 
experimental data. Along these lines, there has been significant progress in recent years 
in perturbative calculations within thermal field theory. The pressure in QCD, for exam- 
ple, is now known to order Unfortunately, an analysis of the convergence of this 
expansion shows that the successive perturbative approximations do not converge for ex- 
perimentally accessible temperatures. This lack of convergence, while not surprising, needs 
to be addressed in order to provide systematic methods for calculating quark-gluon plasma 
observables. 

The lack of convergence of the weak-coupling expansion is not restricted to QCD. In fact, 
even in simple massless scalar field theories similar convergence problems are encountered. 
This indicates that the problem might be universal. The universality of the problem means 
that the technique needed might be quite general and since calculations within scalar theories 
are technically simpler than in full QCD these theories can provide an important testing 
ground for methods to deal with this problem. Like QCD, the weak-coupling expansion for 
the pressure of a massless scalar field theory with a ^^^0^/4! interaction, is known to order 



V = V; 



ideal 



1 - 7« + + ¥ flog + 0.40 

(1) 



4 3 4 V 27rT 



(log ^ - ^ log« - 0.72) a^'^ + Oia' log a 

where Pideai = (7rV90)T^ is the pressure of an ideal gas of free massless bosons, a = 
(y'^(/i)/167r^, and g{jj,) is the MS coupling constant at the renormalization scale fi. In Fig. |l|, 
we show the successive perturbative approximations to P/Pideai as a function of g{2nT). 
Each partial sum is shown as a band obtained by varying fi from ttT to AnT. To express 
g{lj,) in terms of g{2nT), we use the numerical solution to the renormalization group equation 
/Li^Q! = P{a) with a five-loop beta function [^]: 

n a = 3a^ + 32.54a^ - 271.6a^ + 2848.6a^ . (2) 

Ofx 3 

The lack of convergence of the weak-coupling expansion for large coupling is evident in 
Fig. |l|. The band obtained by varying /i by a factor of two is not necessarily a good measure 
of the error, but it is certainly a lower bound on the theoretical error. Another indicator of 
the theoretical error is the deviation between successive approximations. We can infer from 
Fig. |I] that the error grows rapidly when g{27iT) exceeds 1.5. 

A similar behavior can be seen in the weak-coupling expansion for the screening mass, 
which has been calculated to order (7^ . In Fig. 0, we show the screening mass normal- 
ized to the leading order result mLo = 5'(27rT)T/-\/24 as a function of g{2T[T), for each of 
the three successive approximations to m^. The bands correspond to varying yU from ttT to 
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FIG. 1. Weak-coupling expansion to orders 5^, , g^, and g^ for the pressure normalized to 
that of an ideal gas as a function of g{2TrT). 

AttT. The poor convergence is again evident. The pattern is similar to that in Fig. |l|, with 
a large deviation between the order-r/^ and order-gf^ approximations and a large increase in 
the size of the band for g'^. 

There are several ways to reorganize perturbation theory to improve its convergence. One 
method is Fade approximants 0. This method is limited to observables like the pressure, 
for which several terms in the weak-coupling expansion are known. Its application is further 
complicated by the appearance of logarithms of the coupling constant in the coefficients of 
the weak-coupling expansion. However, the greatest problem with Fade approximants is 
that, with no understanding of the analytic behavior of V at strong coupling, it is little 
more than a numerological recipe. 

An alternative with greater physical motivation is a self-consistent approach 0. In this 
method, perturbation theory is reorganized by expressing the free energy as a stationary 
point of a functional fl of the exact self-energy function Il{po, p) called the thermodynamic 
potential 0. Since the exact self-energy is not known, 11 can be regarded as a variational 
function. The "^-derivable" prescription of Baym |^ is to truncate the perturbative expan- 
sion for the thermodynamic potential Q and to determine 11 self-consistently as a stationary 
point of Q. This gives an integral equation for the self-energy that is hard to solve nu- 
merically, unless n is momentum independent. A more tractable approach is to find an 
approximate solution to the integral equations that is accurate only in the weak-coupling 
limit. Such an approach has been applied by Blaizot, lancu, and Rebhan to massless scalar 
field theories and gauge theories [|TOi|TT| . 

Another variational approach is screened perturbation theory (SFT) introduced by 
Karsch, Fatkos and Fetreczky 0]. This approach is less ambitious than the ^-derivable 
approach. Instead of introducing a variational function, it introduces a single variational 
parameter m. This parameter has a simple and obvious physical interpretation as a thermal 
mass. The advantage of screened perturbation theory is that it is straightforward to apply. 
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FIG. 2. Weak-coupling expansion to orders g"^, g^, and g^ for the screening mass normalized to 
the leading-order expression as a function of g{2TrT). 



Higher order corrections are calculable, so one can test whether it improves the convergence 
of the weak-coupling expansion. Karsch, Patkos and Petreczky applied screened perturba- 
tion theory to a massless scalar field theory with a 0^ interaction, computing the two-loop 
pressure and the three- loop pressure in the large- limit. 

In Ref. |T^, a detailed study of screened perturbation theory for a massless scalar field 
theory was presented. The pressure and entropy were calculated to three loops and the 
screening mass to two loops. It was shown that, in contrast to the weak-coupling expansions, 
the SPT-improved approximations converge even for rather large values of the coupling 
constant. In Ref. |13|, the sum-integrals for the three-loop free energy were evaluated exactly 
by replacing the sums by contour integrals, extracting the poles in e, and then reducing 
the momentum integrals to integrals that were at most three-dimensional and could be 
evaluated numerically. The resulting expressions, while truncated in the coupling constant 
were "exact" in the sense that they included contributions from all orders in m/T. In this 
paper we continue the study of screened perturbation theory by performing an analytic 
expansion of the sum-integrals in powers of m/T and demonstrate that the first few terms 
in the expansion give an accurate approximation to the exact SPT result. 

The paper is organized as follows. In section II, we describe the systematics of screened 
perturbation theory. In section III, we calculate the free energy and entropy to three loops, 
and the screening mass to two loops, expanding in powers of m/T. In Section IV, we 
calculate the screening mass to two loops using the m/T expansion. In Section V, we briefly 
discuss the two-loop tadpole gap that generalizes the one-loop gap equation. In Section 
VI, we study the convergence properties of SPT-improved results for the pressure, entropy, 
and screening mass using the m/T expansion. Finally, in section VII, we summarize and 
conclude. Necessary calculational details are collected in four appendices. 
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II. SCREENED PERTURBATION THEORY 



The Lagrangian density for a massless scalar field with a interaction is 



24^ 



+ A£ 



(3) 



where g is the coupling constant and A£ includes counterterms. Renormalizability guaran- 
tees that A£ is of the form 



-AZ9. 
2 ^ 



24 ^ ^ 



(4) 



Screened perturbation theory, which was introduced by Karsch, Patkos and Pe- 
treczky |]12| , is simply a reorganization of the perturbation series for thermal field theory. It 
can be made more systematic by using a framework called "optimized perturbation theory" 



that Chiku and Hatsuda |jT^ have apphed to a spontaneously broken scalar field theory. 
The Lagrangian density is written as 



SPT 



-(' 

2 V 



m 



1^0^ + A/: + A£sPT 



(5) 



Here, £q is the vacuum energy density term, and we have added and subtracted mass terms. 



If we set £q = Q and mj 



we recover the original Lagrangian Eq. 



Screened 



perturbation theory is defined by taking m to be of order unity and mi to of order g , 
expanding systematically in powers of g'^ and setting m\ = m? at the end of the calculation. 
This defines a reorganization of the perturbative series in which the expansion is about the 
free field theory defined by 



free 



-£o + Id, 



1 



2 i2 

m (p 



(6) 



The interacting term is 



2 J 2 
int ^ 2"^!*^ 



A£ 



SPT 



(7) 



Screened perturbation theory generates new ultraviolet divergences, but they can be can- 
celled by the additional counterterm in A£spt- If we use dimensional regularization and 



minimal subtraction, the coefficients of these operators are polynomials in g and 
The additional counterterms required to remove the new divergences are 



m 



-ml). 



AC 



SPT 



-A^n 



Am^ — Ami 



(8) 



Several terms in the power series expansions of the counterterms are known from previous 
calculations at zero temperature. The counterterms Ag"^ and Am^ are known to order P| . 
We will need the coupling constant counterterm only to leading order in a: 



^9' 



2e 



a + ... 



9 



(9) 
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We need the mass counterterms Am^ and Amf to next-to-leading order and leading order 
in a, respectively: 



Ami 



a + 



a + ... 



1 5 

2 



+ ... 



m 



m 



1 • 



(10) 
(11) 



The counterterm for A^o has been calculated to order []T3[. We will need its expansion 
only to second order in a and : 



1 1 

4^ + 8^" + 



5 1 
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1 1 ■ 
4^ + 8^" 
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mim -|- -—mi . 
4e 



m 



(12) 



III. FREE ENERGY TO THREE LOOPS 



In this section, we calculate the m/T expansions of the pressure and entropy density to 
three loops in screened perturbation theory. In performing the truncation m is treated as 
a quantity that is 0{g) and include all terms which contribute to order g^. The diagrams 
for the free energy that are included at this order are those shown in Fig. ^ together with 
diagrams involving counterterms. 




Oa 1a 1b 2a 2b 2c 2d 

FIG. 3. Diagrams for the one-loop (Oa), two-loop (la and lb), and three-loop (2a, 2b, 2c, and 
2d) free energy. 



A. One-loop free energy 

The free energy at leading order in g'^ is 



(13) 



where Aq^^o is the term of order g^ in the vacuum energy counterterm Eq. (|T2D. The expres- 
sion for diagram Oa in Fig. ^ is 



Oa 



(14) 
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with ^ defined in Appendix ^ 

Treating m as 0{g) and including all terms which contribute through 0{g^), we obtain 



(15) 



where Xq and T„ are defined in Appendix ^ and Jq is defined in Appendix 0. 

The resulting expression is logarithmically divergent and the pole in e is cancelled by the 
zeroth order term Aq^^o in Eq. ([l^). The final result for the truncated one- loop free energy 
is 



where rh 



m 
2-kT 



90 

and L = log 



1 - ISm^ + 60m^ + 45m^(L + 7) 



(16) 



B. Two- loop free energy 
The contribution to the free energy of order g"^ is 



(17) 



where Ai£^o and Aim^ are the counterterms of order g"^, respectively. The expressions for 
the diagrams la and lb in Fig. |^ are 



^ la 



The sum-integrals in Eqs. (|T8|) and (|1^) are expanded to the order required: 
1 



(18) 
(19) 



T ^ 2 



Ji + T/i - rn^Xa 



(20) 
(21) 



where J„ is defined in Appendix 0. 

The poles in e in Eqs. (^) and (^Tj) are cancelled by the counterterms in Eq. ([T7|). The 
final result for the truncated two-loop free energy is 



1152 
2A~ 



1 - 12m - I2m^ (L + 7 - 3) + 72(L + 7)771^ 
1 — 6m — Qm^[L + 7) 



(22) 
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C. Three-loop free energy 

The contribution to the free energy of order g"^ is 

= ^2a + -^2b + ^2c + ^2d + ^2^0 + T^^Aam^ + -^^^ [A^m')' 



dm? 



2 (9m^ 



A 2|-^laA 2|*^lbA 2 

Aim H 5-Ai5( H 2"Ai"^i , 



(23) 



where we have included all necessary counterterms. The expressions for the diagrams 2a, 
2b, 2c, and 2d in Fig. R| are 



2a 



16^ 



1 



1 



^2h = --r7;9 



48^ Tpqr P2 + ^2 g2 + ^2 /^2 + ^2 (p + g + p)2 + ^2 
^2c = 4^X^ W+^tg (Q2 + ^2)2 ' 

T -'-4 
^2d = -T"^1 



Expanding in powers of to the appropriate order gives 



^2a = 


-is' 

16 


rj2/2 + J2 J2 + 2T^XJJ2 + T3/2 J2 + 2T/1J1X2 - 2m2rXi J2/2 


^2b = 


48^ 


Xball + T3/ball + 4TJiTsun(0) 






^2c = 


1 2 2 


TJ1J2 + J1J2 + T'/i/2 + TJ2/1 - m2TJ2/2 
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^2d = 


1 4 
-4^1 


X2 + T/2 , 







where /baii, ^Tsun, and Xbaii are defined in Appendices |B|, 0, respectively. 

The poles in e in Eqs. (P^)-(PT|) are cancelled by the counterterms in Eq. 
result for the truncated three-loop free energy is 



(24) 
(25) 
(26) 
(27) 

, (28) 
(29) 
(30) 
(31) 

The final 



9'T' 



T2 



2304(47r)2m 



1 - 2m 



- 12m2 5 + 7L + 37 - 8 log m - 8 log 2 - 4 



C(-i). 



"48(47r)2m 



1 + 2m (L + 7 - 3) - 18m^(L + 7) 



m? 



64m 



1 + 2rh{L + 7) 



(32) 
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D. Pressure to three loops 



The pressure V is given by — JF. The contributions to the pressure of zeroth, first, 
and second order in are given by Eqs. (|16|), (^2]), and (|32D, respectively. Adding these 
contributions and setting £q = {] and m\ = m^, we obtain approximations to the pressure 
in screened perturbation theory which are accurate to 0{g^). 
The one-loop approximation to the pressure is 



ideal 



1 - 15m^ + 60m^ + A5m^{L + 7) 



(33) 



The two-loop approximation to the pressure is obtained by adding Eq. (|2^) with m{ = m . 
Vo+i = Pideaij 1 ~ ^'^ + 15ma + 15m^(L + 7 - 3)a 



-30m^ 



l + 3(L + 7)a 



45m^(L + 7) \ . 



(34) 



The three- loop approximation to the pressure is obtained by adding Eq. (^) with ml = iv?: 



Vo+l+2 = "PidealS 1 - -a + 



59 15^ 5 5C'(-3) ,C'(-1) 

-- + — L + -7 - - J + 5^ ^ ' 



12 4 4' 2C(-3) C(-l) 



a 



15 . 
+ -m 

15 



1- 5 + 37 + 7L-81ogm-81og2-4 



C(-i) 
C(-i) 



a 



a 



-m 



1 -6(L + 7)a 



+ —a 
8m 



(35) 



Note that if we substitute the leading-order result for the screening mass, m 
5'(27rT)T/v^24, in Eq. (p5D, we recover the weak-coupling expansion through order g^. 



E. Entropy to three loops 



Given a diagrammatic expansion for the free energy JF, the entropy density S has a 
diagrammatic expansion defined by 



d 

'S'diag = -^^(T, g, m, mi, /x) 



(36) 



where the partial derivative is taken with all the other variables g, m, mi, and n held fixed. 
The one-, two-, and three-loop approximations to S are then obtained by taking partial 
derivatives of the corresponding expressions for the pressure V. 
The one-loop approximation to the entropy S is obtained by differentiating Eq. (p3D 



where 5ideai = (27rV45)T3. 



'^0 — ^ideal 
2 //I c;^^^3 



- 2 - S 45 ^ 4 

1 —m + 15m m 



15 
Y 



(37) 
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The two-loop approximation to the the entropy S is obtained by differentiating Eq.( 



45 



15 



So+i = cSidcaiS 1 - -a + —ma + —m L + 7 - - a 



4 4 



2 



7 



l + 3(L + 7-l)a 



^45., 
+ -m 



(38) 



The three-loop approximation to the the entropy S is obtained by differentiating Eq. 

5 



o , - / 281 5 15^ C'(-l) 5C'(-3) 



45 . 

H m 

8 

15 
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1 - 
-m 



a 



y + 7L + 37 - 8 log m - 4^^^ - 8 log 2 j a 



l-6(L + 7-l)a 



+ 



25 

32A~ 



-a 



(39) 



IV. SCREENING MASS TO TWO LOOPS 

In this section, we calculate the m/T expansion of the screening mass to two loops. The 
diagrams for the self energy that are included at this order are those shown in Fig. ^ together 
with diagrams involving counterterms. 






1a 2a 2b 2c 

FIG. 4. Diagrams for the one-loop (la and lb) and two-loop (2a, 2b and 2c) self-energy. 

The screening mass is defined by the location of the pole of the static propagator: 

p2 + + n(0, p) = at = -ml , (40) 

where Il{po, p) is the self-energy function. This equation can be solved order-by-order in 
powers of a and mf. The solution at zeroth order in g"^ is simply = m^. 

A. One-loop self-energy 

The self-energy to leading order in g'^ is 

Hi = Hia -ml + Aim^ , (41) 

where Aim^ is the mass counterterm of order a given in Eq. ([T0|). The expression for the 
diagram la in Fig. ^ is 
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nia = Ig' 



1 



2^ Tp P'^ + m'^ 
This diagram is expanded to second order in m^: 

1 



Hi, 



:9 



Ii + Th - m 



(42) 



The pole in Eq. (|4^ ) is cancelled by the counterterm Aim^. The final result for the one-loop 
self-energy is 



24 



Hi 



1 — 6rh — 6rh^ (L + 7) 



— m\ 



(43) 



B. Two- loop self- energy 

The contribution to the self-energy of second order in is 



n2(P) = + n2b(P) + Ii2c + ;f^Aim2 + %Ai/ + Asm^ - Aim^ 



drn^ 



9' 



The expressions for the diagrams 2a, 2b, and 2c in Fig. ^ are 



n2b(P) 



(/22 + rn?)2 ' 
1 1 



Q tqr q"^ + w? + w? {p + Q + ny + m? 



(44) 

(45) 
(46) 
(47) 



The diagrams and II2C sue independent of the momentum P. Expanding to first order 
in m^, we obtain 



n 



2a 



■4^ 



1 



n2c = ^g^ml 



I2 + Th 



(48) 
(49) 



The diagram n2b depends on the external momentum P. The equation ( ^OD for the screening 
mass involves the self-energy at po = 0. To calculate the screening mass to second order in 
g"^, we need the analytic continuation of 11(0, p) to = —rn?. The diagram is calculated in 
Appendix 0. The result is 



n2b(o,p) 

where /gun is defined in Appendix 0. 



Xsun(O) + h 



(50) 
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The poles in Eqs. (^8D-(|50D are cancelled by the counterterms in Eq. (0). The final 
result for the truncated two-loop contribution to the self-energy at = and p 



2 2 

— -m IS 



n2(o,p^ 



l + 2rh 



3L + 7 + 1 - 41ogm - 81og2 - 2 



C(-i) 



+ 



327r^m 



1 2 (L 7) m 



(51) 



C. Screening mass 



Since the dependence of the self-energy on the momentum enters only at order g'^ and 
since the leading-order solution to the screening mass is = m, the solution to the equa- 
tion (1^) to order g^ is simply 



ml = m^ + n(0, 

The result for the one-loop screening mass is 

.0 1 



(52) 



m„ = -a 
' 6 



1 — 6m — 6m (L -|- 7) 



(53) 



The solution to order g'^ is obtained by inserting the sum of Eqs. (^31) and (|5TD into 
Eq (1^). Setting m\ = m?, the result is 



6 



1 
2m 



3L + 7 + l + 41ogm-81og2-2 



C(-i) 
C(-i) 



a — 3m > . 



(54) 



If we substitute the leading-order result for the screening mass, m = gi2TiT)T / y2A, in 
Eq. (^4|), we recover the weak-coupling expansion through order g'^ P]. 



V. GAP EQUATION 

The mass parameter m in screened perturbation theory is completely arbitrary. To 
complete the calculation it is necessary to specify m as a function of g and T. One of the 
complications from the ultraviolet divergences is that the parameters m^, g'^, and m\ all 
become running parameters that depend on a renormalization scale /x. In our prescription 
for recovering the original theory, we must therefore specify the renormalization scale /i* at 
which the Lagrangian (1) reduces to Eq. (|]). The prescription can be written 

^o(/i*) = , (55) 
m2(/ij = m^(/ij = m^(T), (56) 

where m*(T) is some prescribed function of the temperature. 
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(a) [I = 2nT 





1 2 3 4 1 2 3 

g(27iT) g(27iT) 

FIG. 5. Solutions m^{T) to the one-loop tadpole gap equation as a function of g{2'nT) for (a) 



II = 2ttT and (b) /i = m^,. Exact SPT curves are taken from Ref. [13 



The prescription of Karsch, Patkos, and Petreczky for m*(T) is the solution to the one- 
loop gap equation: 



1 



JiK/T)T^- (21og^ + l)m^ 



Their choice for the scale was /i* = T. The function Ji{(3m) is defined as 

e 1 



J^{l3m) = 8/3' 



dk 



In the limit f3m — > 0, this integral reduces to 

MPm) — ^ ^ - 47r/3m - 2 (log ^ - 1 + ^^j {Pm)' 



(57) 



(58) 



(59) 



In the same limit, Eq. (|57| ) reduces to 

. o 1 



= -a 
6 



1 — 6m* — 6m^ (L + 7) 



(60) 



The one-loop gap equation is identical to the one-loop screening mass if we choose m 

rris = m=K 



Various mass prescriptions that generalize Eq. ( p7|) were extensively studied in Ref. |]T3 



In this paper, we confine ourselves to using the tadpole mass which is defined by mf 
I (7^(02). This can also be expressed as a derivative of the free energy 



d 

ml = g'^-^TiT, g, m, mi, /i) 



(61) 



mi=m 



where the partial derivative is taken before setting mi = m. 
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The one-loop expression for the tadpole mass is given differentiating Eq. 

. o 1 



m. = -a 
* 6 



1 — 6m — 6m (L + 7) 



(62) 
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g(2jiT) 



g(2jiT) 



FIG. 6. The one-, two and three-loop SPT improved approximations to the pressure as a 
function of g{2'KT) for ^ = 2'kT. Exact SPT curves are taken from Ref. [13|. 
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I 0.96 




g(27tT) 



g(2jcT) 



g(2jiT) 



FIG. 7. The one-, two and three-loop SPT improved approximations to the pressure as a 



function of g{2TrT) for /i = m*. Exact SPT curves are taken from Ref. |13|. 



In Fig. H, we show the truncated solutions m* (T) to the one-loop tadpole gap equation as 
a function of g{2TTT) for (a) fi = 2'kT and (b) /i = m^,. The solutions have been normalized 
to the leading order screening mass = g{2TiT)T / \/24:. The truncated solutions were 
determined by treating m as a quantity that is 0{g) and truncating at a fixed order in 
g. A g"^ truncation of Eq. (|50D, for example, yields ml = a/6, which corresponds to the 
leading order screening mass. The non-trivial truncations, g^ and g'^, are shown as grey 
dashed lines along with "exact" curves from Ref. |13] which are accurate to all orders in 
m/T. As can been seen from the figure, the gap equation converges very quickly to the exact 
solutions for fi = 2ttT while for /i = they do not seem to be converging. The primary 
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difference between the two scales is that in the case = m* there are additional logm. 
It is possible that these logs need to be further resummed. Note that as the renormalized 
coupling constant becomes larger than g{2nT) ~ 4 the uncertainty due to the variation of 
the renormalization scale /i becomes rather large due to the Landau singlarity present in the 
running of g. For this reason, in all results presented, we restrict ourselves to g{2TiT) < 4. 

VI. SPT-IMPROVED VARIABLES 

In this section, we use the solutions to the tadpole gap equation obtained in Sec. to 
obtain successive approximations to the pressure, screening mass, and entropy in screened 
perturbation theory. 

A. Pressure 

The two-loop SPT-improved approximation to the pressure is obtained by inserting the 
solution to the one- loop gap equation ( |^ into the two-loop pressure (pi]). We can simplify 
the expression by using Eq. (|57|) to eliminate the explicit appearance of logarithms of /i. 
This eliminates all the terms of order a and the expression reduces to 




g(27rT) g(27iT) 

FIG. 8. The one-, two and three-loop SPT improved approximations to the pressure as a 
function of g(27rT) for (a) vrT < ji < AttT and (b) ^m* < fi < 2m*. 

In Figs. 1^ and |^, we show truncations of the one-, two-, and three- loop approximations 
to the pressure for /i = 27rT and fi = m^, respectively. The various truncations are shown as 
grey dashed lines along with "exact" curves from Ref. |]T3[ which are accurate to all orders in 
m/T. As can be seen from the figure, the m/T truncations converge very quickly for the one- 
and two-loop approximations with the final two truncations being virtually indistinguishable 
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from the exact SPT solutions. At three-loops, however, one needs to include all terms up to 
before a reasonable approximation is obtained. We therefore conclude that it is necessary 
to include higher order terms in order to fully converge to the exact SPT result at three- 
loops. Also it appears that the mjT truncations converge better for /i = 27rT, than for 
[I = m*. Despite these caveats, at all loop orders presented here, the highest order m/T 
truncation provide an excellent approximation to the exact SPT results. 

In Fig. ^ we show the one-, two-, and three-loop approximations obtained using our 
truncation to the pressure. The bands shown correspond to the results obtained by 
varying the renormalization scale fi over (a) vrT < fi < AnT and (b) < fi < 2m*. 

This figure demonstrates that the g^ truncations of the pressure yield a convergent series of 
approximations which have very small variations with respect to the renormalization scale. 



B. Screening Mass 

The one-loop SPT-improved approximation to the screening mass is given by the 
solution to the tadpole gap equation (pT|). A two-loop SPT-improved approximation can be 
obtained by inserting the solution to the gap equation for the mass parameter into Eq. (p4D. 

In Fig. ^, we show the g^ truncations of the one- and two-loop approximations to the 
screening mass. The bands shown correspond to the results obtained by varying the renor- 
malization scale /i over (a) tcT < ^ < AnT and (b) < fi < 2m*. One can see from 
this figure that the convergence of the m/T expansion for the screening mass is not as im- 
pressive as for the pressure meaning that higher order truncations are necessary to reliably 
describe the screening mass. Also, we again see that the m/T truncations converge better 
for yU = 2itT, than for fi = m*. 



(a) ttT < ^ < 4nT (b) i m.< M < 2nn, 




g(27iT) g(27iT) 

FIG. 9. The one- and two-loop SPT improved approximations to the screening mass as a 
function of g{2TTT) (a) ttT < f^i < AttT and (b) ^m* < fi < 2m*. 
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C. Entropy 



The one, two-, and three-loop SPT- improved entropies are obtained by replacing m in the 
expressions (^)-(|39|) for Sq, Sq+i, and iSo+1+2 with the solution to the one- loop gap equation. 
As was the case with the two-loop pressure we can use the gap equation to eliminate the 
logarithm L yielding the following expression for the two-loop entropy 



0+1 



5i 



ideal 



15 . 2^ n 45 - 

1 m T + 15m m 

2 4 



(64) 



which is the entropy of an ideal gas of 



This is identical to the one-loop expression Eq. ([ 

particles with mass m. 

In Fig. |T0|, we show the 0{g^) truncations of the one-, two-, and three-loop approx- 
imations to the entropy as a function of g{2TTT). The bands shown correspond to the 
results obtained by varying the renormalization scale n over (a) vrT < fi < AirT and (b) 
|m^, < fj, < 2m^. In both cases the 0{g^) truncation provides an excellent approximation to 
the exact SPT result. Again we see that the m/T truncations converge better for = 27rT, 
than for = m*. 



CO 



T3 



1.00 



^ 0.98 - 



0.96 



0.94 



(a) ttT < n < 47iT 




1.00 



0.98 - 



0.96 



0.94 



(b) im,<n<2m. 




I I Two Loop 
I I Three Loop 







1 



2 

g(27rT) 



2 

g(27rT) 

FIG. 10. The SPT-improved two- and three-loop approximations to the entropy as a function 
of g{2iTT) (a) vrT < fi < AttT and (b) < fi < 2m*. 



VII. CONCLUSIONS 



In this paper, we have continued the systematic study of screened perturbation theory 



from Ref. |]T3[. We applied it to the pressure and the entropy calculated to three loops 
and the screening mass calculated to two loops. By performing an expansion of the sum- 
integrals in powers of m/T, we were able to obtain purely analytical results, without having 
to evaluate integrals numerically. 
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Our calculations show that a truncation of the m/T expansion at 0{g^) is sufficient to 
obtain accurate approximations to the exact one- and two-loop results. For the one- and two- 
loop approximations to the pressure and entropy, the numerical results obtained in and 
the truncated m/T expansions are virtually indistinguishable as can be seen in Figs. |^ and|^. 
At three-loop the 0{g^) truncation provides a reasonable description of the pressure, but it 
seems that higher-order truncations are necessary to provide accurate descriptions. The fact 
that the m/T expansions converge quickly is important since performing the "exact" SPT 
calculations is much more difficult than the m/T expansions. An additional beneffi of the 
m/T expansion method is that the final results can be determined completely analytically. 

In Ref. [0, a generalization of SPT to gauge theories based on hard thermal loop (HTL) 
perturbation theory was proposed. The thermodynamic functions such as the pressure and 
entropy were calculated to one-loop order. The two-loop calculation of the pressure in QCD 
based on HTL perturbation theory requires not only HTL propagators, but HTL vertices as 
well. The exact calculation appears to be very difficult. Expansions like the one presented 
here provide the simplification needed to complete the two-loop HTL calculation |[T^. The 
rapid convergence of the m/T expansion in screened perturbation theory is very encouraging 
in this regard. 
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APPENDIX A: SUM-INTEGRALS 

In the imaginary-time formalism for thermal field theory, the 4-momentum P = (PqjP) 
is Euclidean with = P^+p"^. The Euclidean energy po has discrete values: Pq = 2mTT for 
bosons and Pq = (2n-|-l)7rT for fermions, where n is an integer. Loop diagrams involve sums 
over Pq and integrals over p. With dimensional regularization, the integral is generalized to 
d = 3 — 2e spatial dimensions. We define the dimensionally regularized sum-integral by 




(A.l) 



where 3 — 2e is the dimension of space and /i is an arbitrary momentum scale. The factor 
{e'^ /4:7iY is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, /z coincides with the renormalization scale of the MS renormalization scheme. 
The one-loop sum-integrals that arise in the calculations have the following form 

T 
T 



^logP' 



P (P 



'2\n 



(A.2) 
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Expanding in e to the required order, the specific one-loop sum-integrals needed are 



° " Wt) 45 



1 + 0(e) 



\AtxT) 12 



1+1 2 + 2^-^-^ I e 



.4 C(-i)+ C(-i), 



2e 



(47r)2 UttT/ 



i + 27 + [ ^ - 471 j e + [ y7 + 472 - ^C(3) J + Oie 



(A.3) 
(A.4) 
(A.5) 
(A.6) 



The numbers 71 and 72 are the first and second Stieltjes gamma constants defined by the 
equation 



The specific two-loop sum-integrals needed is 

1 



PQ P^Q\P + QY 



. 



(A.7) 



(A.8) 



It was first calculated by Arnold and Zhai in Ref. The specific three-loop sum-integral 
needed is 



PQR p^qm^p + Q + Ry 



rp4 



24(477)2 V47rT/ 
It was first calculated by Arnold and Zhai in Ref. M. 



(A.9) 



APPENDIX B: INTEGRALS 



We also need some three-dimensional integrals. We choose dimensional regularization to 
regulate infrared and ultraviolet divergences. The integrals are generalized to c? = 3 — 2e 
dimensions of space and /i is an arbitrary momentum scale. 



i 



An 



P 



(27r; 



3-2e 



(B.l) 



The factor {e'^ /Any is introduced so that, after minimal subtraction of the poles in e due to 
ultraviolet divergences, fi coincides with the renormalization scale of the MS renormalization 
scheme. 

The integrals that arise in the calculations have the following form 
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Expanding in e to the required order, the specific one-loop integrals needed are 



Qtt \2mJ 



2e 



'0 

h 
h 



"ivr V2m j 
1 / „ X2. 



1 + 2e 



+ 4 + 0{€^ 



TT 



Svrm \2mJ 
The only two-loop integral needed is 
1 1 1 



1 + -e^ + 0{e') 



qr -|- -|- (p + q + r)^ + 'Ti^ 



(87r)2 \2mJ 



+ 6-81og2 + C(e) 



The only three-loop integral needed is 
1 1 1 



'ball 



pqr -|- + -|- (p + Q + r)^ + 



(47r)3 V2m/ 



6e 



- + 8-41og2 + 0(e) 
e 



(B.2) 

(B.3) 
(B.4) 
(B.5) 



(B.6) 



(B.7) 



APPENDIX C: SETTING SUN DIAGRAM 

The only nontrivial sum-integral required to calculate the self-energy to two loops is the 
sunset diagram, which depends on the external four-momentum P = {po,p)'- 



1 1 
fQR Q2 + R^ + m^{P + Q + + m 

The sum-integral (|C.1|) must be evaluated at po = and p^ = —iv?. The setting-sun sum- 



Xsun(^) i Q2 , ^2 ^2 , ^2 , Q , m2 , ^2 " (^'l) 



integral involves a double sum-integral, so there are three momentum regions. The region 
where both Q and R are hard is denoted by {hh), the region where one momentum is hard 
and the other soft is denoted by {hs), and the regions where both momenta are soft is 
denoted by {ss). The contribution from each of these regions are 



= 0{m) , (C.3) 
I'g! = T^/,,, + 0(m2) . (C.4) 
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APPENDIX D: BASKETBALL SUM-INTEGRAL 



Tbaii is the basketball sum-integral: 

^ JPQR + m? + i?^ + S*^ + ' 

where ^ = -(P + g + i?). 

The basketball sum-integral ( p.l| ) involves a triple sum-integral, so there are 4 momen- 
tum regions: {hhh), (hhs), (hss), and (sss). The contribution from each of these regions to 
order is 



X, 



(?t/is) 



^ = T^4aii + O(m^) . (D.5) 

The {hhh) contribution is given by Eq. ( |A.9| ), while the (hhs) contribution vanishes due to 
Eq. ( |A.8|) . The (sss) contribution is given by Eq. ( p.7|) . 
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